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Abstract 

The shear viscosity n in the van der Waals excluded volume hadron-resonance gas model is consid- 
ered. For the shear viscosity the result of the non-relativistic gas of hard-core particles is extended to 
the mixture of particles with different masses, but equal values of hard-core radius r. The relativistic 
corrections to hadron average momenta in thermal equilibrium are also taken into account. The ratio 
of the viscosity n to the entropy density s is studied. It monotonously decreases along the chemical 
freeze-out line in nucleus-nucleus collisions with increasing collision energy. As a function of hard-core 
radius r, a broad minimum of the ratio n/s ~ 0.3 near r ~ 0.5 fm is found at high collision energies. 
For the charge-neutral system at T = T c = 180 MeV, a minimum of the ratio n/s = 0.24 is reached 
for r = 0.53 fm. To justify a hydrodynamic approach to nucleus-nucleus collisions within the hadron 
phase the restriction from below, r > 0.2 fm, on the hard-core hadron radius should be fulfilled in 
the excluded volume hadron-resonance gas. 

PACS numbers: 24.10.Pa, 24.60.Ky, 25.75.-q 
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I. INTRODUCTION 



One of the important discoveries in Au+Au collisions at BNL RHIC is evidence for strong 
collective effects in the form of transverse flow and asymmetric azimuthal (elliptic) flow. To 
account for this, one needs almost perfect fluid hydrodynamics (i.e. a small ratio, rj/s < 0.2, of 
shear viscosity r\ to entropy density s) starting at the very early stage of nucleus-nucleus (A+A) 



collision (see, e.g., [l, 2, 3|, 4] and reference therein). Perturbative QCD calculations 5] have led 
to a much higher value of the viscosity to entropy density ratio, rj/s > 1 (see Ref. 6[). Based 
on these results, the 'almost perfect fluid' in the deconfmed phase is assumed to be the strongly 
interacting quark-gluon plasma (sQGP) [7]. In Ref. [6j, different possibilities for a structure 
of a minimum of the ratio rj/s at the transition between hadron and quark-gluon phases are 
discussed. An absolute minimum of this ratio was identified with the critical point of QCD 
matter jjj, f|. This conclusion has been based in Ref.0] on an empirical observation in liquid- 
gas transitions for other known substances. In Ref. |8| it was shown that certain field theories 
have the ratio rj/s = 1/47T, and it was conjectured that r]/s > 1/Att for all other substances. 
This restriction appears to be close to the quantum (due to Heisenberg uncertainty principle) 
kinetic lower bound of Ref. [9], rj/s > 1/15, for the gas of quarks and gluons (see Ref. [2]]). 

It is of both practical and fundamental importance to know the ratio rj/s in the hadron phase, 
in particular, near the transition to the sQGP. This will be a subject of the present paper. We 
will use the statistical model with full hadron-resonance spectrum and take into account the 
short-range repulsion effects within the hard-core approximation. Earlier calculations of the 

n 

ratio rj/s were done in Ref. [10J for the pion gas and pion-kaon mixture (see recent pion gas 
results in Refs. ll|, Q, rJ). There were also the calculations of the hadron system viscosity 
within the microscopic transport models: UrQMD in Ref. [yj] and URASiMA event generator 
in Ref. Q. 

During last decades statistical models of hadron-resonance gas (HG) have served as an im- 
portant tool to investigate high energy nuclear collisions. The main obj ect of the study have 
been mean multiplicities of produced hadrons (see, e.g., recent papers [16|, ll?!, ll8|] and references 
therein). An extension to the ideal gas picture based on the van der Waals (VDW) excluded 



volume procedure was suggested in Refs. [19 



20 



21] to phenomeno logically take into account 



repulsive interactions between hadrons. This leads to the well known VDW suppression of 
hadron number densities. However, the particle number ratios are (in Boltzmann approxima- 



2 



tion) independent of the proper volume parameter v, if it is the same for all hadron species. 
Thus, a rescaling of the total volume V leads to exactly the same hadron yields as those in the 



ideal hadron gas. Recently the VDW HG mode 
number fluctuations. It was demonstrated 



has been used in Ref. 



221 ] to calculate particle 



that multiplicity fluctuations are suppressed 
in the VDW HG gas. This suppression is qualitatively different from that of particle yields. 
In contrast to average multiplicities, the suppression of multiplicity fluctuations cannot be re- 
moved by rescaling of the total volume of the system. Thus, the hard-core radius of the hadrons 
can be, in principle, straightforwardly connected with data on multiplicity fluctuations. 

The aim of the present paper is to make calculations of the ratio r]/s in the VDW HG model. 
We present our results for different values of temperature T and baryon chemical potential hb 
along the chemical freeze-out line. In this T — /xg region hydrodynamic expansion within the 
hadron phase is expected. The paper is organized as follows. In Section II we present the 
results for the thermodynamic functions in the multi-component VDW gas. In Section III 
the expression for the viscosity in the VDW HG is obtained. In Section IV the viscosity to 
entropy density ratio is calculated along the chemical freeze-out line. This gives the VDW HG 
predictions for this ratio in central nucleus-nucleus collisions in the hadron phase at different 
collision energies. A summary, presented in Section V, closes the paper. 

II. VDW HADRON GAS 

The VDW excluded volume procedure gives the following transcendental equation for the 



system pressure p [191 ]: 



p = exp ( - ^ ) T (j) , (1) 



T 

where v is the hadron excluded volume parameter, 

v = 4-^rr 3 , (2) 

with r being the particle hard-core radius. In what follows it is assumed to be the same for all 
hadron species. In Eq. (Q]), is the total particle number density in the ideal HG: 

<KT,p B ) = j> = J2 £-2 < TK - ( y) exp (f ) ' (3) 



3 



with 6i in Eq. (}3l) being the ideal gas particle number density of zth particle species. As in 
Ref. 22j, we use the Boltzmann approximation, neglecting small effects of Bose and/or Fermi 
statistics. In Eqs. f|T|3|) . T is the system temperature, ^ the degeneracy factor of ith particle 
species, the particle mass, and fa = bifi B + s^s + Qit^Q the chemical potential due to the ith 
particle charges (bi,Si,qi) — baryon number, strangeness, electric charge. The summation in 
Eq. ([3]) is taken over all hadron species. Finally, K2 in Eq. ([3]) is the modified Hankel function. 
The zth particle number density in the VDW HG equals to: 

n . = exp(- vp/T) <pj _ Xj ^ 
1 + v exp(— vp/T) (ft 1 + v x 

In Eq. (j3J) the notations, Xi = exp(—vp/T)(pi and x = Yli x i> have been introduced. The factor 
R = exp(—vp/T) (1 + v x)' 1 presents the VDW suppression factor. Since the same excluded 
volume parameter for all particle species is assumes, it is the same for all particle densities, and 
the VDW HG energy density. 

The VDW HG entropy density s is equal to: 

dp x f T d(p\ 

1 + vx + — — . (5) 



dT 1 + vx V ch dT 



III. SHEAR VISCOSITY 



The shear viscosity in a non-relativistic gas of hard-core spheres can be calculated analyti- 



cally (see e.g., Ref. [231]): 



5 VmT 

The Eq. ([6]) demonstrates a simple dependence of the viscosity on the system temperature T, 
particle mass m, and hard-core particle radius r. The derivation of Eq. ([6]) from the molecular- 
kinetic theory demonstrates the dependence of the viscosity, 

V oc n I (|p|) , (7) 

on the particle density n, the mean free path I oc l/(nr 2 ), and the particle average thermal 
momentum: 



IpI) 



Jo°° P 2 dp exp[— p 2 /{2mT)] V n 



f™p 2 dpp exp[- p 2 /(2mT)] /8mT 
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The Eq. (JS]) can be easily extended to the relativistic thermal motion: 



„ J °°P 2 dpp exp[- yV + m 2 /T] / 8^T K 5/2 (m/T) 

{lW / °Vdp exp[- v/p 2 + m 2 /T] V it K 2 (m/T) ' U 

At m/T ^> 1, the factor K^/2{m/T) / K 2 {m/T) goes to 1, and Eq. (Q is reduced to Eq. (jSJ). 
In the opposite limit, m/T <C 1, Eq. goes to the well known result for massless particles, 
<|P|> = 3T. 

For the the mixture of particle species with different masses, but with the same hard-core 
radius r, Eq. is transformed into: 

^ 64 r 2 K 2 ( mi /T) n' 1 j 

where rij is given by Eq. (jlj), and n = J2i n i * s the total VDW gas particle number density. 
In Eq. (flQjl . the relativistic correction factors, -Ki = K^^ira^jT^jK^imijT^ from Eq. 
are taken into account. To obtain Eq. (flOl) the same mean free path has been assumed for 
different ith hadron species. This is approximately valid because of the same hard-core radius 
r. The corrections because of different particle masses are neglected. Eq. fflOl) . with the 
summation over all particles and resonance species listed in THERMUS 24|, will be used to 
calculate the hadron viscosity. The ratios rii/n in Eq. (JTD~|) can be substituted by the ideal gas 
values, rii/n = (pi /(f), as we have assumed the same hard-core radius r for all hadrons. The 
VDW suppression factor R is then the same for all hadron species, and it is canceled out in the 
ratios of particle densities. 

Note that the viscosity rj is not proportional to the number of hadron species. If one assumes 
the same hadron masses = m and m 2> T, Eq. fflOl) is reduced back to the case of one particle 
species ()6]). On the other hand, the entropy density s increases with the number of species. 
This fact has been used in recent paper 25] to construct a counterexample to the lower bound 
of rj/s = 1/4-7T assumed in Ref. |8|]. 

In the limit r — > the excluded volume parameter v (j2J) goes to zero. In this limit, the 
thermodynamical functions of the VDW HG, p (pQ), (jlj), and s (jSJ) converge to the corre- 
sponding ideal gas expressions, p ld = T0, n l f = fa, and s ld = + Td(p/dT. On the other hand, 
according to Eq. fflOl) the viscosity rj goes to infinity as rj oc 1/r 2 . Thus, the ratio rj/s diverges 
as r -2 — > oo in the ideal gas limit of r — >• 0. 
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IV. VDW HG MODEL RESULTS 



In this section we present the results of the VDW HG for the ratio r\j s along the chemical 
freeze-out line in central A+A collisions for the whole energy range from SIS to LHC The 
values of T, and 7 S at the chemical freeze-out at different collision energies are presented 
in Table I. They are identical to those values in Table I of Ref. [22 ] . 
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TABLE I: The chemical freeze-out parameters T, fj>B, and 75 in central heavy-ion collisions are 
presented at different cm. energies y / S7v r /v- The ratio rj/s is presented in the VDW HG for three 
values of the hard-core radius, r = 0.1 fm, 0.3 fm, and 0.5 fm. The VDW suppression factor, R = 
exp(—vp/T) (1 + vx) , for these hard-core radiuses is also presented at different freeze-out points. 
The factor K n = K^^im-w /T) / ' K2(m n /T) shows the relativistic enhancement for the pion thermal 
momentum at temperature T according to Eq. ([9]). 



Note that the conditions for average energy per particle, (E) / (N) = 1 GeV [26|, zero value 
of the net total strangeness, S = (this defines fis — ^s(T, ^b)), and the charge to baryon 
ratio, Q/B = 0.4 (this defines jiQ = /ig(T, /i^)), remain the same as in Refs. [27j, l28j. The only 
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tiny difference of T and [ib values from those in Refs. 27J, |28| comes because of the Boltzmann 
statistics approximation used in the present paper and in Ref. The dependence of /j,b on 
the collision energy is parameterized as fl6|. /jb (^/snn) = 1-308 GeV (1 + 0.273 y / Sjv r /v)~ 1 , 
where the cm. nucleon-nucleon collision energy, ^snn, is taken in GeV units. The strangeness 
saturation factor, 75, is parameterized a s \y\ . — 1 — 0.396 exp (— 1.23 T/ jj, B )- Both these 



relations are the same as in Refs. 



22 



27|, 



The resulting ratio rj/s in the VDW HG is 



presented in Fig. 1 as a function of hard-core radius r for several chemical freeze-out points 
from Table I. 
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FIG. 1: The ratio r]/s in the VDW HG as the function of hard-core radius r. The lines correspond to 
different T — fiB freeze-out points from Table I: SIS (^/snn = 2.32 GeV), AGS (^snn = 4.86 GeV), 
SPS (V^w = 17.3 GeV), and RHIC (y/sjfit = 200 GeV). 



The values of ^snn quoted in Table I correspond to the beam energies at SIS (2 A GeV), 
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AGS (11.64 GeV), SPS (204, 304, 404, 804, and 1584 GeV), colliding energies at RHIC 
( v /5^7 = 62.4 GeV, 130 GeV, and 200 GeV), and LHC {^/s^ = 5500 GeV). The excluded 



volume correction factor, R = exp (— vp/T) (1 + vx) , and the entropy density, Eq. ([5 



are calculated using the THERMUS package 24]- The relativistic enhancement factor, K n = 
K 5 /2{m n /T)/ K2(m 7T /T), for average thermal momentum is also presented in Table I at different 
freeze-out points for pions, which are the lightest hadrons. The viscosity is calculated according 
to Eq. ( ITUl) . The Fig. [1] demonstrates a presence of a broad minimum of the ratio rj/s near 
r fa 0.5 fm at high collision energies. For fixed hard-core radius r, the ratio rj/s monotonously 
decreases along the freeze-out line with increasing of collision energy. This is shown in Fig. [2j 
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FIG. 2: The ratio r]/s in the VDW HG along the chemical freeze-out line as the function of ^Jsnn 
at different hard-core radiuses. The solid line presents the results for r=0.5 fm and dashed line for 
r=0.3 fm. 
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The presence of strong collective flow in A+A collisions indicates a rather small viscosity. 
According to Ref. 29[ the following inequality, r]/s < 5 , should be satisfied. As seen from 
Fig. [1], this leads to the restriction, r > 0.2 fm , on the hard-core radius from below in the 
VDW HG. 

The Fig. [3] shows the dependence of the ratio rj/s on the temperature T in the VDW HG 
for charge-neutral system (Q = B = S = and 75 = 1). 
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FIG. 3: The ratio r\js for charge-neutral system {Q = B = S = and 75 = 1) as a function of T. 
We fix the crossover temperature between HG and sQGP as T c = 180 MeV. At T < T c the VDW HD 
results are shown at different hard-core radiuses, the solid line corresponds to r=0.5 fm and dashed 
line to r=0.3 fm. The dotted line corresponds to the lower bound, rj/s = 1/4-7T [8j]. At T > T c , the 
triangle symbol corresponds to Ref. [3l| of the pure gluon plasma in the lattice QCD, and the square 



is the estimate of Ref. 



32( in the quasiparticle model of the sQGP. 
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The VDW HG for r =0.3 fm and r=0.5 fm are presented at T < T c , with T c = 180 MeV 
assumed as the temperature of a crossover between the confined and deconfmed phase. At 
T = T c = 180 MeV a minimum of the ratio rj/s = 0.24 is reached for r = 0.53 fm. Note that at 
high collision energies the baryon chemical potential is small, and a hydrodynamic expansion 
within the hadron phase at T < T c is expected in A+A collisions. 

In the deconfmed phase there were several estimates of the ratio ri/s. The estimates of 

n 



Ref. p based on the perturbative calculations |5j give the large value, rj/s > 1. The non- 
perturbative results correspond to smaller values of the ratio rj/s. In the pure SU(3) gauge 
model on the lattice a shear viscosity to entropy density ratio of rj/s < 1 in the temperature 
region 1.4 < T/T c < 1.8 was found in Ref. [30(] . The same upper bound, rj/s < 1, and the 
estimate rj/s = 0.134 ± 0.033 at T = 1.65 T c were obtained in Ref. [3l|. In the quasiparticle 
description of the sQGP it was found |32| a shear viscosity to entropy density ratio r]/s ~ 0.34 
at T « 1.5T C . 



V. SUMMARY 



In summary, the shear viscosity 77 in the van der Waals excluded volume hadron-resonance 
gas model is calculated. We use the same hard-core radius r for all hadron species. The ideal 
gas limit appears to be a singular one: the ratio of the viscosity rj to the entropy density s 
diverges as r~ 2 at r — > 0. This makes an ideal gas picture inappropriate for any kinetic or 
hydrodynamic calculations. Moreover, the ratio rj/s should be small enough to justify the 
hydrodynamic approach to A+A collisions. Near the chemical freeze-out line this leads to the 
restriction from below, r > 0.2 fm, on the hard-core hadron radius. When the ratio rj/s is 
considered as a function of r, a broad minimum, rj/s ~ 0.3 fm for r ^ 0.5 fm is found for 
the chemical freeze-out at high collision energies. For fixed hard-core radius r, the ratio rj/s 
monotonously decreases along the chemical freeze-out line with increasing collision energy. A 
similar behavior is found for the ratio rj/s as a function of T for the charge- neutral system. 
The ratio r]/s decreases with T, and its minimal value is r]/s = 0.24 at T = T c = 180 MeV 
for r = 0.53 fm. Theoretical estimates of the ratio 77/s in the deconfmed phase require the 
non-perturbative calculations, and they are far from being complete. More efforts are needed 
to clarify the detailed behavior of the ratio rj/s at the transition (crossover) between hadron 
and quark-gluon phases. 
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